DESY 12-225 
DO-TH 12/37 



Light-Cone Distribution Amplitudes of the Ground State Bottom Bciryons in 
HQET 

A. Ali, C. Hambrock, A. Ya. Parkhomenko, Wei Wang 



December 14, 2012 



Abstract We provide the definition of the complete set of light-cone distribution amplitudes (LCD As) for 
the ground state heavy bottom baryons with the spin-parities = 1/2+ and — 3/2+ in the heavy 
quark limit. We present the renormalization effects on the twist-2 light-cone distribution amplitudes and use 
the QCD sum rules to compute the moments of twist-2, twist-3, and twist-4 LCD As. Simple models for the 
heavy baryon distribution amplitudes are analyzed with account of their scale dependence. 



1 Introduction 

Precision tests of the unitarity of the quark mixing matrix remain high on the agenda of flavor physics. They 
allow to pin down the Standard Model (SM) description of CP-violation and may reveal physics beyond the 
SM (BSM). On the experimental side, the two B-meson factories at SLAC and KEK, after approximately 
a decade of their operation, have made a great impact on the origin of CP-violation in the quark sector of 
the SM. Especially in h physics, the i?-factory experiments BABAR and BELLE have concentrated on the 
production and decays of the bb resonances, T{AS) and T(5S), yielding precise measurements of their decay 
products, the B- and the _B-mcsons. Theoretically, inclusive decays of the _B-mcsons, such as the radiative 
and semileptonic decays, are under quantitative control, thanks to the use of the Heavy Quark Effective 
Theory (HQET), allowing an expansion in the inverse 6-quark mass, with perturbative QCD allowing to 
calculate corrections in ag{m}j) to each order in l/mjj. Exclusive semileptonic and radiative decays require, 
on the other hand, a precise knowledge of the non-perturbative quantities - the decay form factors - to be 
analyzed precisely. However, also in such cases, the large energy released in the B — )• hX decays, where h 
is a light meson and X — -f,l'^fr,i^Uf , allows to relate several of these form factors. The resulting large 
energy effective theory (LEET), which has been replaced by a more systematic effective theory, the Soft 
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Fig. 1 The SU{3)p flavor multiplets of the ground-state bottom baryons. The flavor SU {3)p triplet with the spin-parity 
= 1/2+ and the scalar light-quark current (jP = 0"*") is shown to the left, and the sextet with the spin-parity = 1/2"'" 
and the vector light-quark current (jP = 1+) is shown to the right (the similar sextet with = 3/2"'" is not presented 
explicitly). The isospin / and the strangeness S of each baryon state are explicitly specified. 

Collinear Effective Theory (SCET), has been put to good use in these decays, resulting in a number of 
theoretical predictions on different observables in various channels, which are found in global agreement with 
the experimental measurements (see Ref. [I] for a review). 

The attention in flavor physics has so far mostly been focused on the meson sector. Specific processes 
involving bottom baryons, such as rare decays involving flavor-changing neutral current (FCNC) transitions, 
are also potential sources of BSM physics. Experiments at the LHC, in particular the LHC6, are already 
analyzing the copiously produced 6-baryons. As more luminosity is collected, the 6-baryon sector will become 
a wider and quantitative field. Theoretically, these decays are more involved, in particular the non-leptonic 
6-baryon decays, but the semileptonic and radiative 6-baryon decays are tractable. As the starting point 
for the calculations of the transition form factors, a precise knowledge of the quark distributions inside the 
baryons is needed. The LCDAs provide this input in the case of light quarks having high virtuality, which then 
become a crucial input to the calculations of the baryonic transitions in the light cone sum rule approach. 
Baryonic transitions involving a heavy quark have been studied since the beginning of the 90s. However, the 
theoretical precision achieved so far is not comparable to the one attained for the corresponding mesonic 
transitions. So far, most of the heavy baryon distribution amplitudes discussed in the literature [2l|3] are 
the ones which are motivated by various quark models, which are not manifestly consistent with the QCD 
constraints. Hence, it is desirable to make theoretical inroads in this sector, guided by QCD. This effort will 
pay off, as the baryons have the advantage over the i?-mesons of providing access to spin correlations in 
various baryon-to-baryon transitions. A full angular analysis of the semileptonic decays of some 6-baryons, 
such as the decays of the Ah and the f2b, will provide a much more detailed view of the underlying dynamics. 
Of special interest are the semileptonic heavy baryon decays, governed by heavy-to-light form factors. In 
this retrospect, apart from the theoretical analysis based on the heavy quark effective theory [HISIIS], the 
simplification of baryonic form factors in the large recoil limit is exploited in Refs. [HIS] (see Ref. [U1 [TU1[TT| 
for an earlier discussion). In the framework of the soft-coUinear effective theory, it is demonstrated that the 
leading-power heavy-to-light baryonic form factors at large recoil obey the heavy quark and large energy 
symmetries |12j . and in particular the universal function (soft form factor) in Ah ^ A is also calculated 
within the light-cone QCD sum rules in conjunction with the effective field theory [7j. 

In a pioneering paper the complete classification of three-quark light-cone distribution amplitudes 
(LCDAs) of the yl{,-baryon in QCD in the heavy quark limit was given and the scale dependence of the leading- 
twist LCDA was discussed. In addition, simple models of the LCDAs were suggested and their parameters 
were fixed based on estimates of the first few moments by the QCD sum rules method. This analysis can 
be extended to all the ground-state &-baryons with both the spin parities = 1/2+ and = 3/2"'". The 
basic steps and some of the results of such an analysis were presented in a conference proceedings |14j . In 
this paper, we provide the details of the calculations, in particular, concerning the = 3/2"*" 6-baryons. 
Our generalization of the work '13| to the full ground state 5-baryon multiplets provides new input for the 
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eventual calculation of the transition form factors in their decays. In addition, the calculated SU{3) f breaking 
effects presented here are of interest for studying the strange-quark properties in hadronic background in 
general. Especially the S'f,-resonances include three different mass quarks, the bottom, the strange and one 
light u- or d-quark. They provide interesting insights on the nature of the strange-quark propagating in the 
QCD background and give some hints on the validity of the quark condensate models. Comparison with the 
if -mesons and the i^g-meson will show to what extent the condensates are universal quantities. 

This paper is organized as follows: We start by discussing the general properties of the single-heavy 
baryons in the heavy quark limit in Sec. [2] introduce the local interpolating operators of heavy baryons 
in HQET in Sec. [3] and discuss the non-local light-cone operators of these baryons in Sec. |4] Section [5] is 
about the renormalization and the scale-dependence of the various matrix elements introduced earlier. The 
calculations of the correlation functions, which are defined as matrix elements of the non-local and local 
operators are performed subsequently in Sec. |6] followed by the discussion of our non-perturbative model in 
SecjT] Sec. [H] contains our results and numerical analysis, and we conclude with a summary in Sec.[9j 



2 Baryon properties in the HQ limit 

Baryons with one heavy quark Q = c, 6 in HQET are classified according to the angular momentum I and 
parity p of the light quark pair, called diquark. The heavy quarks are non-relativistic particles which decouple 
from the diquark in the leading order of the l/mg expansion. The ground-state baryons = 0) with spin- 
parity are characterized by the spin-parity of the diquark. The spins of the light quarks produce two 
states with jP — 0+ and = l"*". The spin wave- function is antisymmetric for the state with = 0+, 
while Fermi statistics of the baryon state and antisymmetry in color space require an antisymmetric flavor 
wave-function. This results in a baryonic state with isospin 7 = constructed from the light u- and d-quarks 
and the heavy quark Q which is called the ylQ-baryon (the spin-parity is — 1/2+). When the spin-parity 
of the diquark is = 1+, the spin part of the baryon wave- function is symmetric which requires symmetry 
of the wave-function in the flavor space. In the case of light u- and d-quarks beside the heavy quark Q, this 
gives rise to two degenerate states with the isospin 7 = 1, which are called Sq- and Z'g-baryons having 
the spin-parities = 1/2+ and ~ 3/2+, respectively. Inclusion of the s-quark increases the number of 
heavy baryons in the multiplet which (in addition to the isospin) is characterized by the strangeness S. If 
S — —1, there are two baryonic states and S'q with = 1/2+ and H'g-baryon with — 3/2+. For 
S = -2, the baryons with = 1/2+ and = 3/2+ are called Qq and ^2*q. 

The SU{?>)y multiplets of the ground-state bottom baryons are shown in Fig. [T] We denote baryons 
from the triplet and sextet by H^~^ and H^~^, respectively. The corresponding masses of such baryons are 
presented in Tab. [T] and shown in Fig. [2] The important Ahir, Sbir and A^K thresholds are also specified in 
the figure. The spectrum of bottom baryons have been enlarged experimentally, substantially thanks to the 
effort done by the CDF and DO collaborations at the Tevatron collider during the last several years. In Tab.[l] 
experimental measurements are taken from Ref. [15], and we also compare with theoretical predictions (based 
on HQET [T^ and Lattice QCD 17 ) for the masses of the ground-state bottom baryons (in units of MeV). 
Here, A = rn^j — nib and the continuum thresholds sq in HQET for = 4.8 GeV are also given (in units of 
GeV). The baryons lying below the corresponding strong decay thresholds At,TT, Si,tt and Ai,K are of major 
importance, since the rare weak decays can be measured without the immense hadronic background. 

The spinor of the baryon from the = 0+ triplet coincides with the spinor of the heavy quark. In the 
case of the baryons from the = 1+ sextets the situation is more complicated. Both the Rarita-Schwinger 
vector-spinor R^^{v) describing the = 3/2+ baryonic sextet and the spinor of the baryon from the 

= 1/2+ sextet are specific products of the heavy-quark spinor Uj{v) and the polarization vector e'^(w) of 
the light diquark, in which v'^ is the four-velocity of the heavy baryon. This is the effect of spin decoupling 
of light and heavy degrees of freedom in the heavy quark limit. The properties of u-y(v) and e^{v) are shortly 
reviewed below, while detailed information can be found in [18] . Subsequently, the baryonic currents are 
constructed. 
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Isospin 

Fig. 2 (Color online) Masses and thresholds of the b- baryons: SUf{3) triplet consists of A),, HjJ", and (in blue) and the 
others belong to the SUp{3) sextet (in magenta). The important hadronic thresholds Ai,it, Ei,it and Ai,K are indicated by 
horizontal dashed lines. The masses are listed in Tab. [T] HQET calculations have been used in the case of lacking experimental 
input. 

Table 1 Experimental measurements |15) and theoretical predictions (based on HQET |16) and Lattice QCD |17) ) for the 
masses of the ground-state bottom baryons (in units of MeV). Here, A = m j — mi, and the continuum thresholds sq in HQET 

for mi, = 4.8 GeV are also given (in units of GeV). 



Baryon 




Experiment 1151 


HQET [161 


Lattice QCD [T7] 


A 


s 







1/2+ 


5619.4 ± 0.7 


5637+^» 


5641 ± 21+3^^ 


0.8 


1 


2 




1/2+ 


5811.3 ± 1.9 


5809j;7g 


5795 ± letll 


1.0 


1 


3 


^b 


1/2+ 


5815.5 ± 1.8 


5809j;7g 


5795 ± I6I26 


1.0 


1 


3 


^b 


3/2+ 


5832.1 ± 1.9 




5842 ± 26+11 


1.0 


1 


3 




3/2+ 


5835.1 ± 1.9 




5842 ± 26l^g 


1.0 


1 


3 




1/2+ 


5791.1 ± 2.2 


57801^^8 


5781 ± 17^^^ 


1.0 


1 


3 


=■0 


1/2+ 


5787.8 ± 5.1 


57801^1 


5781 ± ntll 


1.0 


1 


3 


^b 


1/2+ 




59031" 


5903 ± 12+11 


1.1 


1 


4 


^f* 
"b 


3/2+ 




5903^79 


5950± 21I2? 


1.1 


1 


4 


a- 


1/2+ 


6071 ± 40 


6036tg} 


6006 ± lOl^g 


1.3 


1 


5 


n* 


3/2+ 




6063l|2 


6044 ± 18l 2? 


1.3 


1 
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The spin and polarization sums are given by 

2 3 
J2u\v)u'{v)=^P+ and "^e*;{v)el{v) ^ -g,,,, + v^v,, (1) 

1=1 i=l 

in which = (1 + f)/2 and the normahzation is given by 

-g^'^ Tr [u{v)e;{v)ei.{v)u{v)] = 6. (2) 

This corresponds to the number of degrees of freedom (3 from the hght quark polarization vector, since it 
has to fulfill the transversality condition Vf^e'^{v) — leaving 3 degrees of freedom, times 2 from the heavy 
quark spinor). The polarization vector is normalized as e^{v)e*^{v) = 1. 
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The heavy quark spin has no fixed direction in the heavy quark limit and can be rotated at wiU. The heavy 
baryon vector-spinor e'^{v)u-y{v) is therefore neither transforming as a spin 3/2 nor as a spin 1/2 quantity, 
since Lorentz invariance is exphcitly broken. One has to restore the proper Lorentz-invariant behavior by 
fixing the heavy quark spinor with respect to the direction of the polarization vector e'^(w). This can be done 
by using the transformation properties of a spin 3/2 Rarita-Schwinger vector-spinor R^^^^j{v) [I8i and [19] : 



(3) 



Following [T5] we define 



(4) 



which fulfill the projection operator conditions = P and P3/2 + P1/2 ~ ^- Now the product of the spinor 
Uj{v) and the polarization vector can be split in proper spin representations 



with 



(5) 



u{v] 



u{v) 



(6) 



A further definition often used is 



V3 



(7) 



with 



(8) 



One can check that B?'\{v) fulfills and R^''^^^{v) is a spin 1/2 spinor, with normalization 



1/2m 



.g^^Tr \r^/^^'{v)R^/^''(v) \ = 4 and - (j^^Tr \R^/'^^'{v)R^/'^'' {v)] = -Tr \r{v)R{v) \ = 2, 



(9) 



corresponding to ([2]), where the 3/2 vector-spinor gets 4 degrees of freedom and the 1/2 spinor gets 2. One 
should keep in mind, that this is a global rotation only to get the proper transformation properties. The 
rotation of the heavy quark spin will not affect the following calculations, which are performed in the heavy 
quark limit. 
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3 Local interpolating operators of heavy baryons in HQET 

The most general three-quark current of the heavy baryon with one heavy quark in the heavy quark limit is 
given by the direct product of the diquark, constructed from two light quarks, and the heavy quark field Q 
due to spin decoupling of the heavy degrees of freedom, discussed in the last section. In this limit, the decay 
constants /^j=o and /]^]=o for the = 0+ triplet of bottom baryons are defined by the matrix elements of 
the local currents PUICT] : 



^abc I 



6"^=(0| [gr(0)C75?^9^(0)] Q^^mnr^iv)) = /^lou.iv). (11) 



The decay constants /^]=i and /^j=i for the jP = 1+ bottom baryons are defined as 



•^"^/nl [qfiO)C^^qUO)] Q^(0)|iJr'(«)) - ^ f^^l.e'^ iv)u^iv), (12) 
[qf{0)Cj'^^q',iO)] Q'^^mHr\v)) = \ /;Le^v)u,iv), (13) 



e 



where a, b, c — 1,2,3 are the color indices, ^ = u.7, 7f = 7^ — '^v^, Q is the effective static field of the 
heavy quark satisfying = Q, the index T indicates a transposition, and C is the charge-conjugation matrix 
with the properties CjJ^C~^ = —7^ and C^5C~^ = 75. In Eqs. (12) and (13l, a factor l/VS is introduced 
so that the decay constants defined here agree with the ones in ^20.,21j . Note that the currents which are 
commonly used in the literature and define the decay constants, for example the two interpolating currents 
for the 17;^ -baryon [H] 

J5 (^) = ^"'^ [qf{x)C^-^ql{x)] (^g,, - I 7*M7t.) Q^i^), (14) 
jg (x) = e'^"^ [qf{x)Cfj'^q',{x)] (^9,, - i 7*^7*.) Q'{x), (15) 



can be easily derived from the currents defining the matrix elements in ( 12 ) and ( 13 ) by applying the 
projectors in Q. 

In QCD sum rule calculations of the LCDAs, the time-ordered product of the non-local and local currents 
(more precisely, the Dirac-conjugated local current) is involved. The most general structure of the Dirac- 
conjugated = 1+ local interpolating current (for the = 0^ current, see [13]) is now given by the following 
linear combination; 

Jr'^ix) = Sabc [9?(x)r'C-iq-f (a;)] g^(x), (16) 

where F' = {At + Bf) 7^ , B — 1 ~ A, and A e [0, 1]. Following the arbitrariness in the choice of the 
local current, i. e., the variation in A, will later be adopted as an error estimate, and we give the results 
for A — 1/2, which corresponds to the constituent quark model current that has maximal overlap with the 
ground-state baryons in the constituent quark model picture, in which all quarks are assumed to be on the 
mass shell |21) . 



4 Non-local light-cone operators of heavy baryons in HQET 

The currents used in the determination of the LCDAs are by definition non-local objects, aligned along 
a light-like direction , with coordinates tiU^ for the i-th quark. For simplicity the baryon is defined in 
the frame of the 6-quark = (5^(0), and the Wilson lines, which are needed to ensure gauge invariance 
of the current, are omitted. Their influence will be discussed whenever necessary. Including the light-cone 
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coordinate dependencies one finds for the non-local light-cone operators that they can be written in the 
heavy quark limit as 







[qiih)Cr''q2{h)]Q^iO) Hl-\v) 



1 



(17) 



where the color indicies are suppressed to simplify the equation as well as the transposition superscript on 
gi(ii). In (17), -T^ symbolizes all possible spinor structures allowed by the Lorentz symmetry and tfj^{ti,t2) 



are the light-cone distribution amplitudes corresponding to the current defined by F^. In total there are 
at maximum eight linear independent structures, n^l, 7'^, n'^^, n''^, —ia^'^Vi,, —ia^'^n^, —inf^a'^Pn^Vp and 
e'"'P''v^npjr'y5 with e°^^^ = 1. 

Since the quarks are aligned along the light-like direction n^, it is convenient for the description of the 
quark dynamics to work in the light-cone basis, where an arbitrary vector is decomposed as 



a'' = - (a_n^ + a+n^) + a'^, (18) 

in which n'' = (1, 0, 0, —1) and — (1, 0, 0, 1) are light-cone vectors with = = 0, n^n'^ = 2, and 
refers to the remaining two space-like dimensions, which are perpendicular to both n'^ and fi^. The scalar 
product of two vectors and 6^ is given by 



The coefficients of the decomposed kinematical vectors 
1 ^n^' 



2 \ v. 
e f n^^ 
2 



(19) 

(20) 
(21) 



have been chosen in a way to fulfill the conditions = 1, {v)£ p,{v) = —1 and the transversality condition 
Vp£^{v) — 0. The light-like vertors are chosen in a way that has no perpendicular components. Accordingly, 
e^(w) and e|'(f) = £^{v) — £'[(«) are called transverse and parallel polarizations, respectively. They fulfill 

the conditions £j|^(u)e^(w) — 0, £\\^{v)e'^^{v) ~ — e^, £j_p(w)e^(w) = — 1 and npe^{v) ~ — eu+ (in the 
following, the reference of the polarization vector and its components to the four- velocity v will be understood 
implicitly). The scalar variable e is the measure for the amount of the diquark parallel polarization. Note 
that the index _L refers to being perpendicular to both vectors n'^ and (and also the four- velocity w^), 
while the index t refers to being perpendicular to only. The vector 



(22) 



is the only possible normalized combination (w^ = —1) of and which fulfills the transversality condition 



0. In this notation, the parallel polarization vector is simply given by eff 



e v^^ 



For comparison, let us recall the set of the non-local currents for the spin-parity — 0+ diquark and the 
corresponding LCDAs given in [T3] with {'0", ip"'", V'^i V""} corresponding to {\P^2, ^ '^i' '^■i} a change 
of the labeling of the decay constants, which are labeled here by the superscript (1) for vector and (2) for 
tensor currents: 



1 



[qiiti)CroMt2)]Qy 

[qi(ti)C^5an„q2it2)] 
\qiih)Cj5q2{t2)]Q-, 



ritut2)flllou^, 



'■=°)=^'(tl,t2)/^lo%, 
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(23) 



where inn 
follows: 



aaun^'^y ■ In a similar way, we introduce the set of the parallel currents and their LCDAs as 



(7i(ii)C^92(i2)] Q7 

qi{tl)C(7nnq2{t2)\ Q-y 



[gi(ii)C92(i2)] i^r') = ^^f(*i'*2)/;,ii 4%, 



v'V 



(24) 



and the set of the transverse currents and LCDAs: 



1 







[gi(ii)C^7^ 92(^2)] 
[qi{h)Cj'[q2ih)]Q^ 
tel(^l)Ccrfiri7± 92(^2)] Qf 

[gi(ti)C^7^ 92(^2)] Q7 



H; 



1 

1 
1 
1 



V'l(^l,t2)/ill 




V'}(ii,i2)/^ii 








V'!(il,t2)41l 





(25) 



The notations 7jJ^ = 7^ ^ 5 (^n^ + ^n'^), {7^,^} — {7^,!^} — have been used. All currents inherit the 
transversality property from the polarization vector. Note, that there is no interpolating current involving 
7y', since 71* = —fiv^. As a result, such currents are linear combinations of the parallel currents defined 
in (24). In a more compact form, the operators can be written as follows: 



qi {tin)Crq2{t2n)Qy 
1 



4^3 



Tr 



r. 



fj^U^ (i^+^T«e^V^l(ii,i2) + ea„s^[l"(<i,t2) + — 'T„,^V'I(ii, ^2) 



+ e/^lw4(ii,i2) + ^/ili757ae"^-"*^^r(ti,t2) 



+ ( ^^^"(ii,i2) + ^^^l(ti,<2) - ^#i'(ti,t2) 



(26) 



in which F represents a Dirac matrix and the decay constants /^j=i and /]^j=i correspond to the vector and 

tensor currents, respectively. Using (26), it is straightforward to calculate the relations of the wave functions 
in any given operator basis. The operators in ^ for the restoration of Lorentz invariance can be applied 
to the currents in (231 — (25) to yield the proper interpolating currents used in, e.g., light-cone sum rule 
calculations of form factors. 

The twist is given by an expansion in w+ (see for comparison |23j). The leading, subleading and subsub- 
leading twists are proportional to 1 and respectively. Similar to the case of light vector mesons 

(see, for example, ^24j) one obtains the twist in the baryon case, listed in Tab. [2] 

The argument for the twist ordering is the following. Based on the choice of basis for v^^ , the orientation of 
the currents with respect to v'^ induces a kinematical difference. The light-cone vectors n'^ and n'^ differ only 
in their direction in the spacial coordinates along the ^-direction. Hence, the characteristic feature is, with 
which three-velocity the baryon moves along the z-axis in the spacial direction of n or n. In conclusion, the 
importance (or twist) can be characterized by V-^. However, this is not necessarily the conventional definition 
of "twist=dimension— spin" . 
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Table 2 Twist of the baryon LCDAs, defined in Eqs. p3l to 1251. 



twist: 


2 


3 


4 


scalar 








parallel 








transverse 









The symmetric LCDAs V'£_L(ti, t2), V'y,j_(ii, ^2) and V'|_L(ii,i2) defined in ([24]) and (|25]) (similar for ([23])) 



are normalized to give the corresponding decay constants /^]=i and /|^]=i in the local limit, in which they are 

satisfying ip^'^"(0, 0) = 1. Performing a Fourier transform {ti — > cui and uji — ulo, UI2 — uuj with u = 1 ~ u) 
leads to the following normalizations: 



C30 00 



duji I (ia;2 V'n''l'"(<^i;'^2) = / ^^d^ / duij}?^'^'^ (uj,u) = 1. 



(27) 











The antisymmetric LCDAs V'y'K^ii ^2) are correspondingly normalized by their first Gegenbauer moments: 



ujduj / du{2u- 1)V'||,1(w,m) = 1, 



(28) 



see Sec. [7] for details of the model functions. Here, LOi are the light-like momenta of the light quarks g^, 
w = wi -|- a;2 is the total light-like momentum of the light diquark system and u and u are the momentum 



fractions of the light quarks. Eqs. (27) and (28) indicate that the LCDAs are constructed in a way, such that 



the energy scale is given by the constants /^=i and /^j=i, while the coordinate dependencies are kept in 



the functions ip. Hence, the two decay constants are determined by the interpolating operators (24) and (25) 
in the local limit, i. e., by the ones given in Eqs. (12) and (13). Splitting the energy dependence from the 



coordinate dependencies leads to the advantage that the decay constants can be calculated at higher orders 
than the more complicated distribution amplitudes. 



5 Renormalization and Scale-Dependence of the Matrix Elements 



A non-relativistic constituent quark picture of the 7?^ -baryon suggests that /f^j* ~ /^j" at low scales of order 



f(i) 



Hi 



'Hi 



1 GeV, and this expectation is supported by numerous QCD sum rule calculations [^Ull^l^l^ . In fact, 
the difference between the two couplings is only obtained at the level of NLO perturbative corrections to the 
sum rules [2HI22j and is numerically small. 

The scale dependence of the couplings is given by 



/l;k.) = /l;k.o)(^ 



asiflo) - asifJ') 7i 



(i) 



47r /3o 

where the anomalous dimensions of the local interpolating operators 
rfln/^](M) 



(29) 



dln/x 



k 



An 



(30) 
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Fig. 3 Types of diagrams for the non-local vertex renormalization. Each of the diagrams a), c), d) and f) appears twice (one 
for each light quark). The thick and thin solid lines correspond to the heavy and light quarks, respectively, and the dashed line 
represents the gauge link (Wilson line). 



are known to NLO [27, and the /3- function expansion 



rfa(/i) 
din /I 



= -/3oa2(^)-/3i a^{p) + 



(31) 



was used to NLO with the coefficients Po — 2 (11 — 2nf /3) and /3i — 4 (51 — 19n//3). 

The decay constants /a, « 0.03 GeV^ and fs, ~ 0.038 GeV^ are known at NLO in the SU{3)f hmit 



while the SU{3)-p breaking effects are only known at LO. The decay constants /f^^jLi and fltj=i coincide 

approximately at the renormalization scale /i ~ 1 GeV ^3j. Note that these couplings cannot coincide at all 
scales since the corresponding operators have different anomalous dimensions. 



f(2) 



The LCDAs are defined via the relations (231 to (25 1. Thus they are not directly measurable quantities 



and obey renormalization group equations. They can be viewed as non-local three-quark vertices, just as the 
local interpolating currents can be viewed as local vertices. The gauge links, also called the Wilson lines, take 
account for the non- locality. Up to a first order expansion of the Wilson lines in gs, the non-local vertices 
are given by: 



fc7 



= -e„-fcr"'3l _ k,+)S{u2 - fc2+), 



^(^2 - k2+) [S{uJi - fci+) - S{uJi - fc4+ - fci+)] 



(32) 



in which color (i, j, ■ ■ ■) and spinor (a, /?,...) indices are explicitly shown together with the light cone momenta 
Wi, ki^ and the heavy quark momentum k. The vertex for the case, in which the gluon line is attached to the 
lower Wilson line is obtained by replacing 1—^2. The expansion of the Green's functions in gg corresponds 
to an expansion of the Wilson lines in gg. In the n-th order expansion term, there are in general n gluon 
fields attached to the vertex, but one-loop calculations are sufficient for the computations in this work. 



and additional vertices to the ones in Eq. ( 32 ) are not necessary. The evolution kernel at one-loop order 



is determined by the ultraviolet poles of the matrix elements of the bare operators. The types of diagrams 
which contribute to the kernel of the renormalization group equation up to one- loop order are given in Fig. |3] 
For the leading twist (neglecting masses) the diagrams give the same contribution for the SU{3)f triplet 
and sextet (Fig. [T]), where only for the transverse distribution amplitude diagram b) vanishes. The SU{3)f 
triplet evolution equation has already been derived in |13j . 

At one-loop order the diagrams shown in Fig. |3] involve at the maximum two quarks in the loop. Hence, 
the derivation of the evolution kernels is almost identical to the meson case. The heavy-light contributions. 
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involving the 6-quark and one light quark, are similar to the heavy-light mesons |26| and the light-light 
contributions are similar to the light-light mesons 27 . The diagram a) is UV finite as it is found in [55]. 
The diagram e) vanishes since the corresponding vertex is proportional to n^n'^ and the gluon propagator 
is proportional to 17^1,. The evolution equation for the leading twist reads il3j: 



dki+ 7^^ (wi , ; /i) i/)" (fci+ , ^2 ; ^) 



+ / rffc2+7 {^^2,k2+;fi)ip'^{uji,k2+;^i)- / dvV{u,v) ip"{vuj,vu}; fi) + -tp"'{uji,uj2; fJ-) 



, (33) 



in which 



(jj 4 



(34) 



is taken from [55] and V{u,v) is the well-known Efremov-Radyushkin-Brodsky-Lepage kernel [551129) : 



V{u,v) 



1 - It 
1-v 



u — V 



0{u 



V — u 



0{v 



(35) 



For transverse leading-twist LCDA -0" diagram b) vanishes and one needs to replace 



V{u, v) -J> V{u, v) 



e{u ~v) + - 0{v - u) 

1 — V V 



(36) 



while the rest stays unchanged. This efi^ect on the evolution is negligible. The "®" and subtractions are 
defined as follows: 



dk+ [7(c., fc+)]^ /(fc+) = J dk+^{uj, k+) [/(fc+) - f{u)] , 


1 

[V{u,v)]^ = V{u,v)~5{u~v) J dtV{t,v). 





(37) 
(38) 



For small evolution steps, ln(^//xo) < 1, the differentiation with respect to /i in Eq. (33) is given by the linear 
approximation: 



d . V (wi,a;2;Ai) - "0 ('^i,W2;mo) 

dfi \n{fi Ho) 



(39) 



and the evolution of the distribution amplitudes can be calculated easily by substituting the initial condition 



■!/'"'(wi, W2; A*o) in the integrand in (33). In the following we give the parallel LCDAs as example. As will be 
shown, the effect of the renormalization is within the errors obtained from the variation of A characterizing 
two different structures in the local interpolating currents. 
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Jr'-,'{x) 







Fig. 4 The correlation function in the QCD background. The light quarks move from the point x in the QCD background 
(shaded regions) to the light-like direction n'^ (dashed line). The b-quark is indicated by the thick black line. 



6 QCD Sum Rules 



Following the standard procedure of QCD sum rules [501151] , we calculate the matrix elements defined in ( 24 ) 



by the local current defined in ( 16 ), i. e 



and (25 1 by approximating the baryonic state 

J(x) |0). The corresponding correlation function is defined by the matrix element 77^/^'^-^' (ti, t2, x) 
i (0| T {Or-y{ti,t2)Jr'-y'(.^)} |0)> which Or'y{ti,t2) can be any non-local operator defined in Eqs. (24) 



and ( |25[ ). The coefficient functions describe the propagation of the quarks inside the baryon from point x to 
the light-cone started from the point 0, which is parametrized by the light quark positions tj, presented in 
Fig. |4] The procedure of constructing the QCD sum rules is well-known and results in the following general 
form: 



3 A- (^/ili 



Rf(2) 



A/t ^ 



= B[nk]{ti,t2;T, So), 



(40) 



where fk — f[l]=i for the even twists and fk = /S=i for the twist-3 parallel distribution amplitudes (and 



swapped according to (25) for the transverse twist). The effective baryon mass is introduced as the differ 
nib, T is the Borel parameter, and sq is the continuum threshold. The r.h.s. in Eq. (40|) 



ence A = m^j=i 



is the Borel-transformed continuum-subtracted invariant function determined through the correlation func- 
tion 7Tfc(ti,t2;r, So). 

The Fourier transform of the correlation function is then given by 



Urr'-fY {^1,(^2] E) = 



dti 
2n 



dt 
2tt 



2 ^iui2t2 



{0\Orj{ti,t2)Jr'Y{x)\0). 



(41) 



Inserting and contracting the quark fields and performing the Fourier transformations yields 



-iEv 



Sm.—r. 



(42) 



The propagators Sq, discussed later, are not free but describe the dynamical evolution of the valence quarks 
within the low-energetic QCD background of the baryon by including condensate contributions. 

The heavy quark spin structure can be seen directly and is given by P+, see Eq. ([T]). The light quark 
polarizations are given implicitly inside the trace by the F and F' structures. The correlation functions must 
be proportional to the Lorentz structure given by the spin sum ([T]). The currents are then automatically 
projected on the proper polarization, such that 



(43) 
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which defines the scalar functions n±{Lu,u] E) and TTy (w, u; i?). They are identical for the parallel and 
transverse polarizations up to different contributions of the local interpolating currents, which gives 

n^{uj, u- E)a,b = i7|| (w, u- E)b,a- (44) 

Since the LCDAs are given for the central value B — 1/2, the distribution functions V-'(w, u : /i) are the 
same for the transverse and parallel polarizations. The dependence of the correlation functions on A and B 



is due to our choice of the set of linear independent operators for the non-local currents in (24) and (25). 
Rewriting the parallel currents by introducing a different basis, for example 7||', would change the above 
described symmetry in A and B. 

The correlation functions are evaluated in the following by using configuration space techniques, instead 
of working in the more common momentum space. The reason is that the leading order sunrise-diagrams 
become very simple in the coordinate space. Sunrise diagrams are certain types of 1 — > 1 multi-loop self 
energy diagrams. (For a review about the evaluation of diagrams of the sunrise type in coordinate space 
see [32]). The non-local diagrams, which are necessary for the calculation of the LCDAs, differ in principle 
from the sunrise type because they do not involve two local interpolating currents, instead they are ripped 
open along a light-like line at the left side of the diagram pictured in Fig. [4j Hence, a more appropriate term 
would be lacerated sunrise diagram. But despite the difference, configuration space techniques work also well 
for this kind of diagrams (as shown below). 

In the following we shortly introduce the different propagators of light and heavy quarks in the configu- 
ration space which are needed for the calculation of the correlation functions. In the heavy quark limit the 
Fourier transform of the heavy quark propagator has the very simple form of a classical point-like particle 
with a non-relativistic on-shell Dirac structure: 

S^=M^^^ [ dt5\x-vt). (45) 



To take the effects of the QCD background for the propagators of the light quark fields into account, the 
method of non-local condensates [331134] is used. The propagator is then given by a sum of the free propagator 
Sq{x) and the universal non-perturbative part Cq{x): 

S,(x) S,(x) C,(x) 



with 



= 2^ " 4^' ^ 12 < (^^^ 

in which color and spin indices are omitted since the color structure is simply given by the Kronecker delta. 
For the non-local condensate we adopt the model proposed in [S5lf5^ : 

CO 

< q{x)q{Q) >=< qq> [ dve'^^" f{v), (48) 







i [a — Z) iiiq 

in which < qq > is the local condensate, rrig is the ratio of the mixed quark-gluon and quark condensate, 
and A is the correlation length. For completeness we also give the Fourier transformed operators {ti — >■ uti, 



uji — u(jj, UJ2 — uuj), which are used in Eq. (42) 



Sq{uuj,t) = -le ^ — ^-^ 



Cq{uuj,t) = e / I • (49) 
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The non-local condensate can be interpreted as the contribution of the quark, which descends to the non- 
perturbative sea of particles populating the QCD background, at one point and reemerges at a different 
point. It is still a question to what extent this sea is universal to the collective of all hadrons. Moreover, little 
is known about the shape of these functions, and different models, such as [S^, have been proposed in the 
literature to describe this non-perturbative behavior. 



Inserting the propagators in the QCD background, Eq. ( 42 ) reads diagrammatically 



n{uj,u;E) 





(50) 



in which the thin lines correspond to the light quarks and the thick line corresponds to the 6-quark. There is 
no 5-quark condensate term, since it is suppressed by l/m^, and hence it vanishes in the heavy quark limit. 
The Borel transform can easily be applied together with the Fourier transform, since for a function f{t) 



iEt 



f{t)dt 



(51) 



holds. As known, the Borel transform handles the UV divergences automatically, making renormalization 
rather easy. As a short reminder, the renormalization condition is applied by adding a certain polynomial 
in E to remove the UV divergent terms, which vanishes due to the derivative of arbitrary order in E in the 
Borel transform. 

The momentum cutoff is applied by cutting off the upper spectrum of the Laplace transform L^^ (mapping 
to the spectral density) of the perturbative contributions, i. e., 



fir) 



e-^/^ fis) ds 



cutoff 



e-^/^ fis) ds, 



(52) 



usmg 



ria) 



dse-'/^ {s-n)"- 



(53) 



One has to be careful, however, what one calls perturbative since there might be mixed perturbative and 
condensate contributions. Mixed contributions originate from diagrams in which one light quark is described 
by the low-energy condensate and the other light quark is described by the high-energy perturbative propa- 
gator. We adopted the procedure of [13] and call the distribution perturbative when at least one perturbative 
term is present. This procedure leads to incomplete /^-functions. For book keeping, we introduce the function 



Ea{x) ^ — ^ fdWe'^ = 1 f , (54) 

r(a + l)Jn r(a+l) 

in which r{a + l,x) — dtf^e"* is the incomplete r'-function. The application of the quark-hadron duality 
thus results in e"^^/'^ r" — > e^^^l'^ t" Ea-\(sQ — ft). For a = iV G M, this function is reduced to the well-known 
form 

^ ri 

En{x)^\-^-'-y^—^. (55) 

^ — ' TV. 

Tl = 

In numerical estimates, it is also helpful to use the relation 

K(a:) = (x) + ^!^^^, (56) 
i (a -I- 2) 
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obtained from (54 1 after integration by parts, connecting the negative value of the parameter a in Ea{x) 
with a corresponding positive one. 

Inserting the propagators (47) and (49) in Eq. (42 1 and performing the Borel and Fourier transforms, the 
sum rules defined in (40) are obtained in a straightforward way. We summarize our sum rule results where 
the leading-twist transverse formula is given as 

\f^. (Af^, + S/^l.) u) e-^'^ = (57) 

3t-4 _^ _ 

[Buj'^ uu + All! {rh2U + rhiu)] Ei{2su:)c " 



r2 



[Aluu + Brh2] /(2™u) E2-a{'is,,) e~ 



3 

[Aluu + Brm] f{2Tuu) E2-a{'^s^) 

/i ~ 

IB 

+ ^ (^i^i) (9292) f{2TLou) J{2tuju) S3-2a(2sKs) e-", 

in which = sq - a;/2, m = 77i/(2t), uj = i^/(2t), = S(^/(2t), = - k/2, = s^^ - k/2, 
Skr = Suj — k/2 — k/2, and the short-hand notations 



A A 



(58) 



/S- [Af^^.+Bf^^^,)^l{u:,u)e-^'^ (59) 
{[Aw + B{rhi+ m.2)] E2{2s^) + Buj {m2U + rhiu) Ei{2s^)}e~'^ 



2UUJT ' 2UUJT 

are used. For the twist-3 and twist-4 LCDAs, the results can be derived analogously: 

lf(l) ( AfW ^uA^) "iJ,! 

_ 

~ 47r4 



[B ^3-a(2sK) + (Bciu + 2A7JI2) S2-a(2s,)] /(2™m) e- 
('Z2g2>r3 + (^Boju + 2Ami) E2-ai2s-^)] fi^ruju) e" 



27r2 
2A 

+ (9191) (9292) /(2rwM) /(2™u) £;3_2a(2s,B) e-"", 



3 

~ 47r4 
S(gi9i>T3 



{[AcZ; (u - u) + B (mi - 7712)] i?2(2S(^) + (to2U - TOi?2)i?i(2s^)} e"" 

[^3_,(2S«) - WU^2-a(2sK)] /(2TWU) C"" 



27r2 

[^3-a(2ss)-L:;^/i?2-a(2ss)]/(2™M)e-'^, 



^/i'r (^/iJi- +^410 (61) 
[B ^3(2s^) + A (mi + m2) £;2(2s„)] e"" 



3 

~ 27r4 



[v4 S3_„(2s,) + Sm2 ^2-a(2s«)] /(2™u) e-"^ 
^^"^ [Ai?3-a(2ss) + Smi£;2-a(2ss)]/(2™u)e-'^ 
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Fig. 5 Model functions for the fe-baryon LCDAs, composed of the exponential part for the heavy-light interaction and the 
Gegenbauer polynomials for the light-light interaction. 



Table 3 Moments for the LCDA model functions in terms of their parameters. 



twist 


(1> 


{^-') 










2 


ao 


ao/3eo 


3ai/5 


ai/5ei 


802/7 


02/7^2 


twist 


(1> 












3 


ao 


ao/2eo 


ai 


ai/2ei 




a2/2£2 


4 


ao 


ao/eo 


a I 


ai/ei 


12 


12 /e2 



2S 

+ (9191) (9252) /(2rwM) /(2™u) £;3_2a(2s«s) e"". 
The result for the parallel counterpart can be obtained by the replacement A-<^ B. 



The QCD sum rules given in (57) can not be used in calculations directly. The main reason for this is that 
the sum rules are built from a patchwork of different contributions, the perturbative and condensate parts. 
They show neither smooth behavior, nor necessarily the correct asymptotic behavior, i. e., the asymptotic 
behavior of the perturbative contribution. As a consequence, one has to propose model functions which are 
then constrained by the sum rules. The LCDA models are discussed below in Sec. [Tj 



7 LCDA Models 

In the parametrization, in which the baryon is described by the total energy of the light quark system uj and 
the momentum fractions u and w, the following arguments become manifest. The dynamics of the heavy- 
light baryon part, i. e., the dynamics of the heavy quark and the light diquark, is described with the same 
formalism which has worked well for the heavy-light meson dynamics. The dynamics of the light quarks in 
the diquark is then described in the same way as the light-light mesons. This ansatz is consistent with the 
equations of motion for the light quarks |24| . In conclusion we choose the multiplicative ansatz from |13| . In 
this approach the Gegenbauer polynomials {x) |38j take the light quarks into account and an exponential 
factor characterizes the dynamics of the heavy-light system, as shown in Fig. [5] It is more convenient to 
choose Cn^'^ix) Gegenbauer polynomials for twist-2 but C}/'^{x) for the other twist functions similar to the 
expansion of the vector mesons 24J . The first three polynomials are sufficient to account for the precision in 
this work. They are defined as follows: 



C^(x) = 1, C^{x)^2\x, C^{x) = 2\{l + \)x'^ - \. 



(62) 
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Fig. 6 The dependence of the transverse twist-2 model parameters listed in Tab.jsjon the variation of the coefficient A of the 
local interpolating currents (the parallel ones are obtained by the replacement A — > 1 — A). 

Table 4 Input values from Refs. |13II31II39II40| . For A and sq, see Tab. Ill 



X 





128 MeV 
0.16 GeV^ 
0.6 GeV 



m 

(ss)/(qq) 



-242 MeV^ 
0.8 
0.8 GeV^ 
1.36 GeV2 



To obtain the model fit, we calculate the momentum fraction integrals (the moments) which are defined for 
an arbitrary function f{uj,u) as: 



2S() 



The model functions for the LCD As of different twists are: 



a-n CfJ^{2u - 1) 



n=0 



\a 



3/2, 



2 

n—0 



1/2, 



7 . s a« Cy^(2u - 1) 



Tl = 



\a 



1/2, 



in which the twist is indicated by the subscript numbers and 

r-l 







u - 



1) 



with C, 



1/2 I 2 



3/2 I 2 



1, \Ci 



1/2|2 



^3/2|2 



3, C. 



1/2 I 2 



1/5, and [Cj 



3/2 I 2 



(63) 

(64) 
(65) 
(66) 

(67) 

6. The pref actors in 



1/3, \Ct 

front of the sums in Eq. (64) — (66) {lu'^u{1 — u), etc.) are determined by the corresponding perturbative 
part in order to give the correct asymptotic behavior (compare the B term in the first line in Eq. (57)). The 
parameters are strictly positive to satisfy the asymptotic behavior. The moments of the functions (64 1 — 
( 66 ) which are calculated with the use of ( 63 ) are hsted in Tab.jsj 

In the construction of the models for the LCDAs, we have truncated the Gegenbauer expansion at the 
second non-asymptotic term and have taken the limit sq ^ oo in the integral over uj which has negligible 
effect. 



8 Results and Numerical Analysis 



The moments of the functions defined in Eqs. (64) — (66) can be calculated using (63). To perform a 
numerical analysis, we discuss and specify the required input parameters which are summarized in Tab. [4j 
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Table 5 Parameters for the model functions in Eqs. | |64[ l — ( |66| l for the transverse LCDAs (the parallel ones are obtained by 
the replacement A I — A and the Ai, and Si, currents correspond to the parallel Sf, and ones, respectively). For twist-2 
LCDAs, the parameter dependence on A is plotted in Fig. [6] The twist notations 3a and 3s correspond to the LCDAs ip"" and 
i/i-", respectively. 



a2 



eo[GeV] £i[GcV] 



£2[GeV] 



2 

3s 
3a 
4 



6.4A 1.4A+0.6 
A+0A4 A+5.7 
0.12A-0.08 0.56A-0.77 



A-IA 



A-2.6 



0.35A-0.43 



-0.07A-0.05 0.65A+0.22 
A+0.34 A+1 



A-0.17 
0.25A-0.16 
A+0.41 



5.5A+3.8 
A+29 



twist 



<ll 



0.2 



El 



£2 



2 

3s 
3a 
4 



0.25A+().46 
A+0.68 

0.04A-0.14 
A-1.6 



-0.16A+0.16 



6.6A+0.6 
A+0.68 
0.12A-0.09 

A-1.6 
0.17A-0.17 



A+6.7 
0.56A-0.91 

A-2.9 
O.llA-0.11 



0.03A+0.11 
A+0.16 



-O.lA-0.03 0.63A+0.38 



0.57A+1.1 

A+4 
-27A+92 

160 
0.39A-0.49 

A-1.3 
-0.82A-3.1 



A+0.61 



A+1. 3 



0.36A+0.03 

A-0.02 
0.3A-0.24 

A+0.54 
0.33A-0.33 

A-1 
1.2A+0.34 
A+4.1 



twist 



ao 



a2 



£0 



El 



£2 



2 

3s 
3a 
4 



8A+1 1.3A+1.3 
A+1 A+6.9 

0.17A-0.16 0.56A-1.1 
A-2 A-3.22 



-O.lOA-0.01 0.62A+0.62 
A+1 A+1. 62 



0.45A-0.63 



0.41A+0.06 
A+0.11 

0.44A-0.43 
A+0.27 



0.87A+0.07 
A+2.53 



The values of the effective baryon masses A = m^j — mi, in the HQET for mt — 4.8 GeV are presented 

b f— — ■ 
where experimental measurements |15j and theoretical predictions (based on the HQET |16j and 

Lattice QCD [17 ) for the masses (in units of MeV) of the ground-state bottom baryons are also shown. A 

comparative analysis of predictions for the heavy-baryon masses can be found in Refs. |16|l41j. The continuum 

threshold values Sq (the last column in Tab.[l]) used by us are in agreement with the ones from [TB] derived for 

the baryon-mass evaluated to order l/nit within the HQET. For the discussion of these parameters see [5^ 



and references therein. Note that the shape function f(i/) in the non-local quark condensate (48) is assumed 
to be flavor independent for all light quarks. 

The calculation is performed in the 6-baryon rest frame with u+ = 1 at an energy scale of = 1 GeV. 
The method of the non-local condensates which involves the parameters A and uiq is not yet completely 
understood. Especially, since there is only one model parameter known, namely the ratio 

of the 5-dimensional and 3-dimensional local condensates. This parameter determines the center of the quark 
virtuality distribution in the QCD background, but is not sufficient to determine the shape of the quark 
distribution. To determine the form, also yet unknown dimension-7 local condensates are needed. We took 
the shape parameter A as the universal parameter which is not influenced by either the baryon or the mass 
of the propagating quark. For the strange quark, where 

2 (sD^s) 

"4 = (69) 

the situation is more difficult since even dimension-5 condensates are not yet clearly understood. The value R, 
which is defined by 



(70) 
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Fig. 7 (Color online) The Twist 2 functions of E^, (blue, solid), Hj, (red, dashed) and flf, (cyan, dot-dashed) at the energy 
scale /i = 1 GeV (blue, solid line) and the energy scale ii = 2.5 GeV (blue, dashed line) including the most conservative error 
A G [0, 1] for X'i, (blue, shaded region). 




Fig. 8 Model function ip(uj,u) of the 6-baryon H;,. Shown are the twist 2 (upper left hand), symmetric twist 3 (upper right 
hand), antisymmetric twist 3 (lower left hand) and twist 4 (lower right hand) LCDAs. All entires are evaluated at the energy 
scale ^ = 1 GeV. 



varies from R sa 0.8 [TOH^ which gives mg ~ m^g to values around R w 1.3 [i51HilH5] which gives 

1.7 ml, (71) 



which we took for our calculation since in this case the 5'/7(3)F-breaking effects appear already in the lower 
w-spectrum. The influence of the choice of R on our results is within the already large uncertainties due to 
the other parameters. With Tab. [3] the model parameters are extracted from the moments and collected 
in Tab. [sjin which the first moment is 1 by definition according to ([27]), and the are required to be non- 
negative. The twist 2 model parameter dependence on A is given in Fig. [6] as example. To give an overview 
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over the tables of model parameters, we show the plots of the LCDAs in Fig [7] in the w and u distributions. 
The corresponding LCDAs for Eb (S'^) in the (w, u) plane are shown in Fig. p^as example. 

9 Conclusions 

As shown in Fig. [T] we are able to obtain the LCDAs for the entire ground-state multiplets of the bottom 
baryons with one heavy quark, thereby generalizing the work by Ball, Braun and Gardi |13j . We accounted 
for the mass-breaking effects due to the strange-quark mass and calculated the evolution of the LCDAs with 
the use of the renormalization group equations, obtained from the one-loop renormalization of the non-local 
light-cone operators given in Eqs. (24) and (25 1. The resulting LCDAs are plotted in Figs, [t] and |8] in which 



the evolution from /i = 1 GeV to fi — 2.5 GeV and our error estimates are also shown. We find that the 
5[/(3)F-breaking effects are of order 10 percent. The sources for the SU{3)f breaking are the strange-quark 
mass, the Borel parameter r, the momentum cutoff sq and the non-local condensates. The latter ones are 
not well determined in the case of the strange quark. However, the mass-breaking effects appear in the 
lower part of the w-spectrum in the region of the constituent s-quark mass, as expected. Even though the 
S'[/(3)F-breaking effects appear to be quite large at some momenta, they are within the conservative errors 
for the massless case, as is the evolution to the energy scale /i = 2.5 GeV. The errors are obtained by varying 
the parameter A in the full possible range [0,1], which defines the linear superposition of the two local 



interpolating currents, given in Eq. (16 1. Reducing the variation of A to the less conservative range [0.3, 0.7] 



results in a variation of the distribution amplitudes of roughly the same size as the evolution to /i = 2.5 GeV. 
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